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Abstract 

In his 1685 paper "Observationes cyclometricae" published in Acta Eruditorum, Adam 
Adamandy Kochariski presented an approximate ruler-and-compass construction for rectifi- 
cation of the circle. It is not generally known that the first part of this paper included an 
interesting sequence of rational approximations of 7t. Kochariski gave only a partial explana- 
tion of the algorithm used to produce these approximations, while promising to publish details 
at a later time, which has never happened. We reconstruct the complete algorithm and discuss 
some of its properties. We also argue that Kochariski was very close to discovery of continued 
fractions and convergents of n. 

1. Introduction 

Adam Adamandy Kochariski SJ (1631-1700) was a Polish Jesuit mathematician, inventor, and 
polymath. His interest were very diverse, including problems of geometry, mechanics, and astron- 
omy, design and construction of mechanical clocks, perpetuum mobile and mechanical computers, 
as well as many other topics. He published relatively little, and most of his mathematical works 
appeared in Acta Eruditorum between 1682 and 1696. He left a reach correspondence, however, 
which currently consists of 163 surviving letters yj. These letters include correspondence with 
Gottfried Leibniz, Athanasius Kircher SJ, Johannes Hevelius, Gottfried Kirch, and many other 
luminaries of the 17-th century, giving a rich record of Kochariski's activities and a vivid descrip- 
tion of the intellectual life of the period. Recently published comprehensive monograph [2] gives 
detailed account of his life and work, and includes extensive bibliography of the relevant literature. 

Among his mathematical works, the most interesting and well-known is his paper on the recti- 
fication of the circle and approximations of K, published in 1685 in Acto Eruditorum under the title 
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Figure 1: Kochariski's construction of approximation of n. 

Observationes cyclometricae ad facilitandam praxin accomodatae Annotated English trans- 
lation of Observationes with parallel Latin version has been made available online by the author 

a. 

The paper has three distinct parts, the first one giving a sequence of rational approximations of 
K. This will be the main subject of this note, so more about this will follow in the next section. 

The second part of Observationes is the one which is most often commented and quoted. There 
the author proposes an approximate solution of the problem of the rectification of the circle, giving 
an elegant and simple construction of a linear segment whose length approximates K. Figure 1 
shows this construction exacly as Kochahski had it the original paper. We start with drawing a 
semi-circle of radius AB centered at A, inscribed in a rectangle BGHD. Than we draw a line AI 
such that ZMC = 60°. When the lower side of the rectangle is extended so that HL is equal to the 
diameter of the circle, and a line is drawn from / to L, one can easily show that 



|/L| = ^y^l20- 18^3 = 3.1415333..., (1) 

which agrees with K in the first four digits after the decimal point. This compass and ruler con- 
struction is often referred to as Kochanski's construction. 

The third part of [3| gives yet another approximation of n, this time expressing it as a sum of 
multiples and fractional parts of 1/32, 

96 4 1 1 1 3217 

+ JL + i . _L + _^ = ffil = 3. 1416015625. 
32 32 2 32 32-32 1024 

It is fair to say that if the rectification of the circle reported in Observationes received a lot of 
attention from both contemporaries of Kochahski and historians of mathematics [I5l[6ll3, then the 
first part of his paper has mostly been forgotten. In what follows we will show that this is perhaps 
unjustly so, as it includes some intriguing sequences of fractions approximating n, origin of which 
has not been explained by commentators of Kochariski's work. 
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2. Sequence of rational approximations of Ti 



In the table on p. 395 of [3J (also p. 2 of [4J), Kochahski gives the following sequence of pairs of 
lower and upper rational approximants of it: 



25 221 r333 3551 f 1667438 16677931 f 9252915567 92545833601 
T'T j ' 1 106 ' 113 j ' t 530762 ' 530875 j ' \ 2945294501 ' 2945825376 J 

136727214560643 1367364691440031 
43521624105025 ' 43524569930401 J ' 



(2) 



As mentioned in the original paper, two of these fractions can be further reduced, 5307^52^ = 25538 1 ' 
2945825376 ~ 30685681 ' '^^ile all the Others are already written in their lowest terms. He then 
partially describes the algorithm generating these approximants, which could be explained using 
modern terminology and notation as follows. Let us denote the first element of the n-th pair (lower 
approximant) by P„/ Q„, and the second element (upper approximant) by Rn/Sn- The approximants 
are then generated by recurrence equations 

Qn+l=SnXn + l-, (3) 
Pn+l=RnXn + ^. (4) 
Sn+l=Sn{Xn-^l) + \, (5) 

i?„+l =i?„(x„ + l)+3, (6) 

where Rq = 22, 5o = 7. In these formulae Xn is a sequence of numbers which Kochahski calls 
genitores, giving the first four values of x„: 

15,4697,5548,14774. (7) 

Unfortunately, he does not explain how he obtained these numbers. He only makes the following 
remark regarding them: 

Methodicam praedictorum Numerorum Synthesin in Cogitatis, & Inventis Polymath- 
ematicis, quae, si DEUS vitam prorogaverit, utilitati publicae destinavi, plenius ex- 
ponamj^ 

In spite of this declaration, Cogitata & Inventa Polymathematica have never appeared in print. It is 
possible that some explanation could have been found in unpublished manuscripts of Kochahski, 
but unfortunately all his personal papers gathered by the National Library in Warsaw perished 
during the Warsaw Uprising in 1944. To the knowledge of the author, nobody has ever attempted 
to find the algorithm for generating the sequence of genitores. In the most comprehensive analysis 
of Kochahski 's mathematical works published up to date (H, Z. Pawlikowska did not offer any 
explanation either. 

In the subsequent section, we attempt to reproduce the most likely method by which Kochahski 
could have obtained the sequence of genitores, and, consequently, a sequence of approximants of 
K converging to n. We will also explain why he gave only first four terms of the sequence. 



will explain the method of generating the aforementioned numbers more completely in Cogitata & Inventa 
Polymathematica, which work, if God prolongs my life, I have decided to put out for public benefit, (transl. H.F.) 
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3. Construction of genitores 



It seems plausible that the starting point for Kochahski's considerations was Archimedes approxi- 
mation of ;r by 22/7 and the result of Metius 

333 355 

It is also likely that Kochahski then noticed that Metius' result can be obtained from Archimedes' 
approximation by writing 

333 _ 22-15 + 3 

106 ~ 7-15 + 1 ■ ^ ^ 

Where is the factor 15 coming from? The key observation here is that 15 is the "optimal" factor, 
in the sense that it is the largest integer value of x for which ^^^^^ remains smaller than n. 

The next most likely step in Kochahski's reasoning was the observation that the upper approx- 
imant can be obtained by incrementing 15 to 16, 

355 ^22^ (10) 
113 7-16 + 1 

Q C C 

Repeating this procedure for ^ produces another pair of approximants, 

355-4697 + 3 _ 1667438 

113-4697 + 1 ~ 530762 ' 
355-4698 + 3 _ 1667793 

113-4698 + 1 ~ 530875 ' ^ ^ 

where 4697 is again the largest integer x for which ^^^'^^^ < n. Recursive application of the above 
process produces the desired sequence of pairs given in eq. (|2]), and the values of x thus obtained 
are precisely what Kochahski calls genitores. 

What remains to be done is proving that the above algorithm indeed produces a sequence of 
lower and upper approximants of n, and that these converge to K in the limit of n — )■ oo. 

4. Kochanski approximants 

We will present the problem in a general setting. In what follows, a will denote a positive irrational 
number which we want to approximate by rational fractions. 

Suppose that we have a pair of positive integers R and S such that their ratio is close to a but 
exceeds a,R/S > a. Together with [aj , we then have two rational bounds on a, 

ifi<a<| (13) 

Suppose now that we want to improve these bounds. As we will shortly see, this can be 
achieved by considering fractions which have the form 
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where x is some positive integer. Before we go on, let us first note that the function f{x) = , 
treated as a function of real x has positive derivative everywhere except at x = —1/5, where it is 
undefined, and that there exists x where f{x) = a, given hy x= {a — [a\ )/{R — aS) . 

Definition 1 Let a be a positive irrational number, and let R and S be positive integers such that 
I > oc. Genitor ofR,S with respect to a will be defined as 



ga{R,S) 



a - [a\ 
R-aS 



(15) 



Let us note that if ga{R, S) is positive, then it is the largest positive integer x such that < a, 

i.e., 

ga{R.S) = max{x e N : -^^^ < «}■ (16) 

In this notation, the four genitores given in the paper can thus be written as gjt{22,l) = 15, 
^;r(355, 113) = 4697, 1667793, 530875) = 5548, and g;r(9254583360, 2945825376) = 14774. 
Using the concept of genitores, we can now tighten the bounds given in eq. ( [T3| ). 

Proposition 1 For any a G IQ^ and R,S G if ^ > (X and if the genitor gaiR^S) is positive, 
then 

^ Rga{R,S)+[a\ R{ga{R,S) + l) + [a\ ^R 

Sga{R,S) + l Siga{R,S) + + 1 



The second and third inequality is a simple consequence of the definition of gaiR^S) and eq. ( 16). 
The first one can be demonstrated as follows. Since R/S > a, then R/S > [a\, and therefore 
Rga{R:S)>[a\Sga{R:S). Now 

Rga{R,S)+la\ > [a\Sga{R,S) +la\, 

and 

Rga{R,S) + [a\ 
Sga{R,S) + l >L«J' 
as required. To show the last inequality let us note that 

R{ga{R,S) + \)+[a\ R [a\S-R 



S{ga{R:S) + 1) + 1 S S{S{ga{R:S) + 1) + 1) ■ 
Since R/S > [a\ , the numerator is negative, and the last inequality of ([Tt]) follows. □ 



The above proposition gives us a method to tighten the bounds of (13), and the next logical 
step is to apply this proposition recursively. 

Definition 2 Let a ElQ^ and let Rq^Sq be positive integers such that Ro/So> oc and gaiRQiSo) > 
0. Kochariski approximants of a starting from Rq^Sq are sequences of rational numbers {P„/ 2„},7=i 
and {Rn/Sn}'^^Q defined recursively for n G N U {0} by 

P„+i=RnXn+[a\, (18) 

Qn+ 1 — SfiXn + 1 , 

Rn+1 =R„iXn+l)+[a\, 

Sn+l = (JC„ + 1 ) + 1 , 
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where x„ = ga{.RniSn)- Elements of the sequence {Pn/ Qn\'^=i will be called lower approximants, 
and element of the sequence {i?„/S,j}^^Q - upper approximants. 

Note that 

P„=R„-R„i, (19) 

Qn^Sn-Sn^U (20) 

therefore it is sufficient to consider sequences of R„ and 5„ only, as these two sequences uniquely 
define both upper and lower approximants. 

Proposition 2 Kochanski approximants have the following properties: 

(i) Xn is non-decreasing sequence of positive numbers, 

P] R Rq 

(ii) [a\ < < a < < — for all n > I, 

R,j 

( Hi) — is decreasing, 

(iv) is increasing, 

(v) lim — = lim — = a. 



For (i), because of the definition of x„ = gaiPmSn) shown in eq. ( 15), we need to demonstrate 
that Rn — OcSn is non-increasing. To do this, let us check the sign of 

Rn - aSn - {Rn+\ " C(S„+l) = Rn - aSn - {Rn{Xn + 1) + ["J ) 

+7t{Sn{x„ + \) + l) = a-[a\-{Rn- aS„)xn 

a - [aj 



a-[a\- {Rn - aSn) 



Rn CiSn 



The last expression, by the definition of the floor operator, must be non-negative, thus i?„ — nSn 
is non-increasing, and Xn is non-decreasing as a result. Now, since the definition of Kochanski 
approximants requires that xq is positive, all other Xn must be positive too. 

Property (ii) is just a consequence of the Proposition 1, which becomes clear once we note that 

Rn ^ i?„-i(x„_i + l) + 3 

Sn Sn-\{Xn-\ + l) + l' 

and 

Pn _ Rn-lXn-l +3 
Qn Sn—lXn—l + 1 

where Xn- 1 = ga {Rn- uSn-i). 
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To show (iii), let us compute the difference between two consecutive terms of the sequence 

Rn_ _ Rn \ _ Rn-\yn-l + [O^J _ Rn \ 

where we defined _y„_i = ga{Rn-iiSn-i) + 1- This yields 

Rn_Rn-l^ (Rn-iyn-l + L«J)-S'n-l - R„-l{S„^iyn-l + 1) 
1 Sfi—i 
_ [Ct\S„^i — Rn-\ 



Sn-l{Sn-iyn-l + 1) 



<o, 



because, by (ii), Rn-i/Sn-i > [(x\ . The sequence Rn/Sn is thus decreasing. Proof of (iv) is similar 
and will not be presented here. 

Le us now note that Rn/Sn is bounded from below by a and decreasing, thus it must have a 
limit. Similarly, ^ = ^" g"'' is bounded from above by a and increasing, so again it must have 

a limit. To demonstrate (v), it is therefore sufficient to show that limits of and J' are the 
same, or, what is equivalent, that 

We start by defining jn = Rn/Sn — Pn/Qn and observing that 

Rn Rn~ Rn—l Rn—l^n Rn^n—l 



rn = 



Sn Sfi Sfi—l Sn{Sn Sn—l) 



By substituting Rn — Rn-i{xn-i + 1) + [ocj and Sn — 5'„_i(jc„_i + 1) + 1, one obtains after simpli- 
fication 

Rn-,-[a\Sn^, _ |fr-L«J 
Sn{Sn-Sn-l) 5„(/--l)- ^ ^ 

Since Rn/Sn is decreasing, and starts from Rq/Sq, we can write 

, , §-^«J |-^«J (23) 

where we used the fact that Sn = Sn-i{xn-i + 1) + 1 and where = gaiRn-\,Sn-i)- Since Xn 
is non-decreasing, and Sn increases with n, we conclude that 7„ — )■ ad n — oo, as required. □ 

Let us remark here that x„ is indeed only non-decreasing, and it is possible for two consecutive 
values of x„ to be the same. For example, for a = a/2 and Rq /Sq = 3 /2, we obtain 

{xX=o = 2, 4,4, 15, 17,77, 101, 119,..., 

where xi=X2. 
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5. Initial values 



One last thing to explain is the choice of the starting values Rq, Sq. Definition [2] requires that the 
genitor of these initial values is positive, so how can we choose Rq.Sq to ensure this? We start by 
noticing that the second pair of Kochahski's approximants (Metius' fractions y^, yj|) are known 
to appear in the sequence of convergents of the continuous fraction representation of Tt. As we 
shall see, this is not just a coincidence. 

Let us first recall two basic properties of continuous fraction expansion of a positive irrational 
number a, 

1 

a = aQ-\ . 

ai H 



1 

a2H 



1 

«3H 

a4 + ■ • . 

By convergent Pn/^n we will mean a fraction (written in its lowest terms) obtained by truncation of 
the above infinite continued fraction after a„. The first property we need is the recursive algorithm 
for generating convergents and values of a„. 

Propositions Consecutive convergents Pn/Qn of a can be obtained by applying the recursive 
formula 



OCqn-l — Pn-l 



(24) 



Pn - OCqn 

Pn+l= Pnan+l+Pn-U (25) 
<?«+ 1 = 1 + <?n- 1 , (26) 

with initial conditions ao = [a\, ai = [ ^i^^ J, pq = ao, qo = \, pi= aoai + 1, = ai. 
For example, for a = ;r we obtain 

3 22 333 355 103993 104348 208341 312689 



\Pn 


r -\ 


[qn , 





1' 7 ' 106' 113' 33102 ' 33215 ' 66317 ' 99532 ' 

833719 1146408 4272943 



265381' 364913 ' 1360120' 



(27) 



Convergents are know to be the best rational approximations of irrational numbers, which can 
formally be stated as follows. 

Proposition 4 Ifpn/qn is a convergent for an irrational number a and p/qis an arbitrary fraction 
with q < qn+i, then 

\qna-pn\ < \qa-p\ (28) 

Elementary proofs of both of the above propositions can be found in ||9l. We also need to recall that 
convergents Pn/qn are alternatively above and below a, so that for odd n we always have p„/qn > 
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a, and for even n, Pn/qn < oc. Suppose that we now take some odd convergent pik+i/ qik+i, and 



further set p = [a\ , q = I. Inequality (28 1 then becomes 

P2k+i - aqik+i <a-[a\, (29) 

and hence 

a — I a I 

^ ^ > 1. (30) 



Pik+i — ccqik+i 



This, by the definition of the genitor given in eq. ( 15 1, yields gaiPik+i^qik+i) > 0. leading to the 
following corollary. 

Corollary 1 Ifp2k+i/ qik+i is an odd convergent of a, then ga {pik+i, qik+i) > 0, and Rq = pik+x, 
Sq = q2k+\ can be used as initial values in the construction of Kochanski's approximants. In 
particular, one can generate Kochanski's approximants starting from the first convergent of a, by 
taking Rq — a^ai + \, Sq = a\, where ao = [a\, ai = [l/(a — ao)\. 

Note that Kochariski in his paper indeed started from the first convergent of n, by taking Rq/Sq = 
22/7. Obviously, if one starts from the first convergent Rq = pi, Sq = qi, then the first lower 
approximant will be the second convergent. Pi = p2, Q\ = qi, and indeed in Kochanski's case 
P\/ Q\ = Pi/ q2 = 333 / 106. Other approximants do not have to be convergents, and they normally 
aren't, although convergents may occasionally appear in the sequence of lower or upper approx- 
imants. For example, in the case of a = ;r, = p-i/q^ = 355/113 and P2/Q2 = Ps/qs = 
833719/265381. 

We should also add here that the choice of the first convergent as the starting point is the most 
natural one. Among all pairs Rq^Sq where Sq < 106, the only cases for which gj:{Ro,So) > are 
R = 22k, R = Ik, where k E {1,2, ...,15}. If one wants to obtain fractions expressed by as small 
integers as possible, then taking = 1 is an obvious choice. 



6. Concluding remarks 

We have reconstructed the algorithm for construction of rationals approximating n used in [3|, and 
we have demonstrated that it can be generalized to produce approximants of arbitrary irrational 
number a. Under a suitable choice of initial values, approximants converge to a. 

Using these results, we can generate more terms of the sequence of genitores for a = n, 
Rq/Sq = 22/1 , going beyond first four terms found in Kochanski's paper: 

{xnTn=^ = {g^{Rn,Sn)}:^Q = {15,4697,5548, 14774,33696,61072, 111231, 

115985,173819,563316,606004,...}. (31) 

We propose to call this sequence Kochanski sequence. It has been submitted to the Online Ency- 
clopedia of Integer Sequences as A191642 filOil . and its entry in the Encyclopedia includes Maple 
code for generating its consecutive terms. 

Knowing that jc„ = ^Z^^g , we can also understand why only four terms of the sequence 
are given in the paper. In order to compute x„, one needs to know n with sufficient accuracy. For 
example, 20 digits after the decimal point are needed in order to compute xq to X3 . Kochanski was 
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familiar with the work of Ludolph van Ceulen, who computed 35 digits of n, and this was more 
than enough to compute X4. Nevertheless, Kochahski in his paper performed all computations 
keeping track of "only" 25 digits, and this was falling just one digit short of the precision needed 
to compute X4. 

It is also interesting to notice that the recurrence equations in Definition |2] strongly resemble 
recurrence equations for convergents Pn/qn in Proposition [3] Kochahski was always adding 3 and 
1 to the numerator and denominator in his approximants, because, as remarked earlier, he noticed 
that 

22-15 + 3 _ 333 22-16 + 3 _ 355 

7-15 + 1 ~T06' 7-16+1 ~ m' ^ 
He apparently failed to notice that 

333-1+22 355 



106-1+7 113' 



(33) 



that is, instead of finding the largest x for which {22x + 3)/ {Ix + 1) < ;r, one can take the last two 
approximants, 22/7 and 333/106, and then find the largest ;c such that (333x + 22) /(106x + 7) > Tt. 
If he had done this he would have discovered convergents and continued fractions. His genitores 
would then be a„ values in the continued fraction expansion of Tt. In the meanwhile, continued 
fractions and convergents had to wait until 1695 when John Wallis laid the groundwork for their 
theory in his book Opera Mathematica fr\\. 

One little puzzling detail remains, however. If we look at the Definition [2} we notice that the 
sequence of lower approximants Pn/Qn starts from n = I, not from n = 0, as is the case for the 
upper approximants Rn/Sn- Indeed, Pq.Qo are not needed to start the recursion. Nevertheless, in 
the table of approximants given in [j3|, in the second row there is a pair of values corresponding 
to n = 0, namely Pq/Qq = 25/8 (in the first row of the table he also gives the obvious bounds 
3 < ;r < 4). These numbers are not needed in any subsequent calculation, and Kochahski does 
not explain where do they come from. One can only speculate that perhaps he wanted the table 
to appear "symmetric", thus he entered some arbitrary fraction approximating n from below as 
Po/Qo- 
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